CELLULAR NETWORK HANDOFF MODELING 



Cross-reference to Related Applications 

[0001] The present application claims the benefit of provisional patent application: Serial 
No. 60/413,449 to Mark et al., filed on September 26, 2002, entitled "Cellular Network 
Handoffs," which is hereby incorporated by reference. 

Background of the Invention 

[0002] The present invention relates to the field of cellular network systems. In particular, the 
present invention provides for efficient modeling and analysis of cellular handoff decision 
mechanisms used in cellular networks. 

[0003] In cellular networks, each mobile station (MS) maintains connectivity via an active 
set of base stations (BS). A handoff algorithm determines the dynamics of the active set as 
the MS moves through the network. In hard handoff, the MS is "handed off* from one BS to 
another BS as it leaves the cell coverage area of the first BS and enters that of the second BS. 
In this case, the active set of an MS consists of at most one BS at any given time. Hard 
handoff algorithms are used in the GSM (Global System for Mobile Communications) and 
GPRS wireless networking standards and are still under active investigation for use in High 
Data Rate (HDR) services. 

[0004] Wireless technologies based on CDMA generally employ soft handoff, whereby the 
MS maintains an active set that may contain multiple BSs. A soft handoff occurs whenever a 
BS enters or leaves the active set of an MS. Soft handoff algorithms are used in the IS-95 and 
cdma2000 standards and have been proposed for use in the WCDMA standard. 
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[0005] The handoff behavior of mobile units in a cellular network may be critical to the 
overall performance of the network in terms of quality-of-service, resource utilization, and 
signaling load. Ideally, the network should maintain seamless quality-of-service for an active 
mobile user engaged in a call as it traverses cell boundaries. Even when the mobile station is 
in the idle state, the choice of cell assignment set may impact the network in terms of resource 
utilization and signaling load, as well as the quality-of-service experienced by the MS when it 
transits to the active state. Thus, the proper design and dimensioning of the handoff algorithm 
may be crucial to the deployment of a cellular network. 

[0006] What is needed is a discrete-time framework for analyzing the performance of 
handoff algorithms executed by an MS based on pilot signal strength measurements from 
candidate BSs. 

[0007] Vijayan and Holtzman were among the first to propose an analytical model for 
handoff based on signal strength measurements. (See R. Vijayan and J. M. Holtzman, "A 
Model for Analyzing Handoff Algorithms," IEEE Trans, on Vehicular Technology, vol. 42, 
pp. 351-356, Aug. 1993.) They modeled the trajectory of a mobile station as a continuous- 
time process wherein handoffs corresponded to the level-crossing events of the relative signal 
strength between two candidate base stations. They further argued that the level-crossing 
events could be modeled approximately by Poisson processes with time-varying rate 
functions. A limitation of the Poisson level-crossing model is that it is accurate only for 
relatively high level-crossings. Furthermore, when the model is applied to the realistic 
scenario of sampled signal measurements, the sampling interval must be sufficiently small for 
the model assumptions to hold. Another limitation of this model is that it is applicable only to 
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relative signal strength measurements and may not accommodate constraints on absolute 
signal strength. 

[0008] Subsequently, Zhang and Holtzman proposed an alternative approach to analyzing 
handoff based on the Gaussian properties of the sampled and processed received signal 
strengths. {See N. Zhang and J. M. Holtzman, "Analysis of Handoff Algorithms Using Both 
Absolute and Relative Measurements," IEEE Trans, on Vehicular Technology, vol. 45, pp. 
174-179, Feb. 1996.) Their approach incorporated absolute signal strength thresholds in the 
handoff algorithm. Approximations for the handoff probabilities are obtained by making 
some simplifying assumptions. However, one skilled in the art will observe that the 
approximate formulas disclosed in this article may be inaccurate, particularly when relative 
signal strength alone is used as the basis for handoff. What is needed is a handoff model that 
can utilize an exact discrete-time model based on the Gaussian properties of received signal 
strengths. 

[0009] In many papers related to handoff analysis, the MS trajectory under consideration is 
limited to the straight line directly connecting the two candidate base stations involved in the 
handoff. However, handoff behavior may depend strongly on the trajectory followed by the 
MS. For arbitrary MS trajectories, simulation-based methods have been the only alternative 
for obtaining accurate performance measures. What is also needed is a system that can utilize 
a representation of a general trajectory in a cellular network capable of yielding a concise 
characterization of handoff performance over a wide range of MS trajectories in the network. 

[0010] Early work on handoff analysis has largely been based on computer simulation 
studies. In industrial practice, computer simulation remains the primary means for choosing 
key parameters to optimize the performance of modern-day wireless networks. Detailed 
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computer simulations of wireless cellular networks require considerable computation time, 
making them cumbersome to use for the purposes of network design and dimensioning. 
[0011] What is needed is a system capable of modeling the efficiency of handoff 
mechanisms in various cellular networks. 

Brief Summary of the Invention 

[0012] One advantage of the present invention is that it provides a discrete-time formulation 
for handoff problem(s). 

[0013] Another advantage of this invention is that it characterizes the handoff behavior of a 
given mobile station in terms of cell assignment events associated with a discrete-time 
random process determined by the averaged measurements of signal strengths received from 
base stations. 

[0014] Another advantage of this invention is that it provides an efficient numerical 

procedure to compute the handoff performance measures. 
[0015] Yet a further advantage of this invention is that it uses an exact discrete-time model 

based on the Gaussian properties of received signal strengths. 
[0016] To achieve the foregoing and other advantages, in accordance with the invention as 

embodied and broadly described herein, the present invention includes a discrete-time 

framework for analyzing the performance of handoff algorithms executed by an MS based on 

pilot signal strength measurements from candidate BSs. 
[0017] In yet a further aspect of the invention, the handoff mechanisms may be primarily 

hard handoff mechanisms, which may be applicable to networks based on the GSM and 

GPRS standards. 
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[0018] In yet a further aspect of the invention, the handoff mechanisms may be primarily soft 
handoff mechanisms, which may be applicable to networking standards such as IS-95 and 
cdma2000. 

[0019] In yet a further aspect of the invention, the handoff mechanisms may be both hard and 
soft handoff mechanisms. 

[0020] Additional objects, advantages and novel features of the invention will be set forth in 
part in the description which follows and in part will become apparent to those skilled in the 
art upon examination of the following or may be learned by practice of the invention. The 
objects and advantages of the invention may be realized and attained by means of the 
instrumentalities and combinations particularly pointed out in the appended claims. 

Brief Description Of The Several Views Of The Drawings 
[0021] The accompanying drawings, which are incorporated in and form a part of the 

specification, illustrate an embodiment of the present invention and, together with the 

description, serve to explain the principles of the invention. 
[0022] Figure 1 is a diagram showing a straight-line trajectory between two base stations in a 

cellular network. 

[0023] Figure 2 is a diagram showing assignment regions based on relative signal strength. 
[0024] Figure 3 is a diagram showing assignment regions based on relative signal strengths 

with absolute thresholds. 
[0025] Figure 4 is a diagram showing nonlinear assignment regions. 
[0026] Figure 5 is a diagram showing handoff as a crossing level problem. 
[0027] Figure 6 is a diagram showing a mobile trajectory in a cellular network. 
[0028] Figure 7 is a diagram showing a reduced geometric structure. 
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[0029] Figure 8 is a diagram showing a minimum geometric structure. 
[0030] Figure 9 is a diagram showing straight-line trajectories within a minimum geometric 
structure. 

[0031] Figure 10 shows probabilities of being assigned to station 55, along the trajectory Trl 
withh= 1. 

[0032] Figure 1 1 shows probabilities of being assigned to base station 55, along the 

trajectory Tr2 with h = 1 . 
[0033] Figure 12 shows probabilities of being assigned to station 55, along the trajectory Trl 

with h = 3. 

[0034] Figure 1 3 shows probabilities of being assigned to station 55, along the trajectory Tr2 
with h= 10. 

[0035] Figure 14 shows probabilities of handoff from 55, to BSj along the trajectory Trl with 
h = 3. 

[0036] Figure 1 5 shows probabilities of handoff from 55, to BSj along the trajectory Tr2 with 
h = 3. 

[0037] Figure 1 6 shows probabilities of handoff from 55, to 55, along the trajectory Tr2 with 
h= 10. 

[0038] Figure 17 shows analysis and simulation for the mean number of handoff following 
the trajectory Trl. 

[0039] Figure 1 8 shows analysis and simulation for the mean number of handoff following 
the trajectory Tr2. 

[0040] Figure 1 9 shows delays in handing off following the trajectory Trl . 
[0041] Figure 20 shows delays in handing off following the trajectory Tr2. 
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[00421 


Figure 21 shows 


crossover noints vs mean number of handoffs for Trl. 
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Fieure 22 shows 

X IgUl V ±* Jllv TT kj 


crossover noints vs mean number of handoffs for Tr2. 
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Fienre 21 shows 


mean number of handoffs h = 1 
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Ficnire 24 shows 


mean number of handoffs h = 3 
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Figure 25 shows 


a cellular network scenario for soft handoff. 


[0047] 


Figure 26 shows 


a grid structure for soft handoff evaluation. 


[0048] 


Figure 27 shows 


a flow diagram of an embodiment of an aspect of the present 



invention. 

[0049] Figure 28 shows a block diagram of an embodiment of an aspect of the present 
invention. 

Detailed Description of the Invention 

[0050] The handoff mechanism employed in a cellular network may have a significant 
impact on overall network performance. Evaluation of handoff mechanism performance has 
for the most part been carried out using brute-force computer simulation or simplified 
approximations. The present invention utilizes a new discrete-time framework for evaluating 
the performance of handoff mechanisms such as hard handoff algorithms based on pilot signal 
strength measurements. Analytical expressions and associated numerical procedures for 
evaluating the cell assignment and handoff probabilities of a mobile station moving along a 
linear trajectory are disclosed. By approximating an arbitrary trajectory with a piecewise 
linear path, a complete characterization of the overall handoff performance of a cellular 
network under a given set of parameters may be obtained. This invention provides an 
accurate and efficient means for the design and dimensioning of high performance handoff 
algorithms. 



[0051 ] The starting point of the present invention is a discrete-time formulation of the 
handoff problem. The formulation characterizes the handoff behavior of a given mobile 
station in terms of cell assignment events associated with a discrete-time random process 
determined by the averaged measurements of signal strengths received from base stations. 
From the formulation, expressions may be derived for the cell assignment and handoff 
probabilities at each sampling epoch along a trajectory traversed by the MS. In the present 
best mode of the invention, the trajectory is an arbitrary straight-line trajectory. However, one 
skilled in the art will recognize that other trajectories may be utilized. A numerical procedure 
may then be developed to compute handoff performance measures. 
[0052J To characterize handoff performance for arbitrary trajectories in a cellular network, a 
general path may be approximated by a piecewise linear path within a reduced geometric 
structure derived from the cellular network geometry. A concise characterization of the mean 
number of handoffs over a wide range of mobile trajectories in the network may then be 
obtained. This characterization preferably provides a measure of the overall signaling load 
incurred by the handoff algorithm. 
[0053] A basic signal propagation system model for a general class of handoff algorithms 
will be discussed. This discussion will first focus on a familiar continuous-time setting and 
then on a discrete-time setting, providing the basis for subsequent discussion. The signal 
propagation model is primarily for the pilot signal strength used by handoff algorithms. The 
propagation model preferably provides the basis for analyzing the handoff performance of an 
MS traveling along a straight-line trajectory at constant speed in a cellular network. A 
discrete-time version of this model may be used in subsequent analysis of handoff 
performance. 
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[0054] The cellular network consists of a set of base stations. The i th base station, BS h 
located by a position vector b iy lies approximately at the center of its associated cell i . The 
coverage area of cell i may be determined by the pilot signal strength from BSi. Assuming 
that the MS is initially located at position a at time 0 and connected to BSu the MS may reach 
position b at time 7, moving at a constant speed of v along the line segment between a and b. 
As the MS moves along this trajectory, the candidate BS for handoff as shown in Figure 1 is 
BSj. The MS may make handoff decisions based on measurements of the pilot signal 
strengths received from the two base stations BSi and BSj . The position, r(t), of the mobile 
user at time t may be given as follows: 

r(t) = a + vte\ 0<t<T 9 

where e is the unit vector in the direction of the vector b - a. 

[0055] The pilot signal strength (in dB) received from a BS as a function of distance in 
typical urban mobile environments may be modeled empirically as the sum of three 
components: path loss, shadow fading, and fast fading. The path loss may be a deterministic 
function of the distance d from the BS, while the shadow and fast fading components may be 
modeled as random processes. This model has been used in a number of studies of handoff 
performance in cellular networks. Under the assumption of constant MS speed along the 
straight-line trajectory, the pilot signal strength maybe equivalently expressed as a function 
of a continuous-time parameter t. 

[0056] We will now describe a continuous-time model. The pilot signal strength (in dB), 
received at the mobile from BSi at time t may be represented as follows: 

Y i (t) = m i (t) + W^t) + Z i (t)u>Q, (1) 

where 



m ) .(0 = v I .- Ml .log(||r(0-6 i W||) > ( 2 ) 
and [W t (t)},i € C are independent, zero mean stationary Gaussian processes. The term rm(t) 
accounts for path loss, while W® models the effect of lognormal (shadow) fading and Zjt) 
models fast fading. The autocorrelation function of Wjt) has been experimentally observed to 
have the form: 



Rw i (T) = a 2 w i exp\ 



■o J 



(3) 



where aw, is the standard deviation of the shadowing signal strength and the constant d 0 is 

called the decay factor. 
[0057] To eliminate the effect of fast fading, the MS applies an exponential averaging 
window to the measured signal strength Yi(t). The impulse response of the averaging window 
may be given by 



1 f ^ 

/"(')=-r exp 



vt 



,t>0, (4) 



d avJ 

where d av determines the effective size of the averaging window. The processed pilot signal 
strength from base station BS t may then given by 

where * denotes the standard convolution operation. 
[0058] Handoff performance based on a discrete-time model for the signal strength 

measurements will now be discussed. In the discrete-time model, the mobile unit may sample 
the pilot signal strengths at time instants t k = kt s , where t s is the sampling interval. The 
distance between successive sampling positions of the mobile may be denoted by d s = vt s . The 
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discrete-time model may capture the handoff behavior of an MS more accurately than its 
continuous-time counterpart since all signal strength measurements are sampled. 
[0059] In discrete time, the basic signals of the system model may simply be sampled 
versions of their continuous-time counterparts: 

Y\k] = Y^m\k] = m^ (6) 

[0060] The processes { W s [k]} ,i € C , are independent, zero-mean, stationary Gaussian 
processes characterized by an autocorrelation function which may be given by: 



Rw i (m) = a 2 w i cxp 



V d o J 



(7) 



One skilled in the art will recognize that the cross-correlation between received signal 
strengths may easily be incorporated into this model. 
[0061] The shadow fading process {WifkJ} cyan be represented as a first order 
autoregressive (AR) process by the following difference equation: 

WifkJ = aWifk - 1] + VJkJ; k>0, (8) 
where (Vifkjj is a zero mean, stationary white Gaussian noise process with variance a* . The 

parameters, (a,cr v 2 ) , of the AR model preferably determine an autocorrelation function for 

{WifkJ} of the form 



[0062] By comparing (9) and (7), the AR parameters for shadow fading may be determined 
as follows: 



a = exp 



V d *J 



(10) 
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< =(!-«>)< , (ID 
[0063] The discrete-time equivalent of the exponential smoothing window may be given by 



r -kd ^ 



, Jfc>0 (12) 



[0064] The processed pilot signal strength from base station BSi may then be obtained as 
follows: 

*,M=/-v[*HM> 03) 

[0065] where * is used here to denote discrete-time convolution. Note that the measured 
signal strength Y t [k] maybe expressed as: 

Y i [k]=m i [k] + W i [k] (14) 
[0066] The relative processed signal strength between base stations BSi and BSj may be given 

by 

X[k]=Xi[k]-Xj[k]. (15) 
[0067] The evolution of the process {X£k]} may be governed by the following second-order 
difference equation: 

X\k] = {a^b)X\k-\}-abX\k-2} + j-(m\k]-am\k-^ k>2 (16) 



av 



x, M = j- {i bm i [°] + m > M) + ( a+b ) v > M + v > W) ' (1 7) 

x i [o]=^(H[o]+^[o]). (18) 

From (16), we see that /XA/j is a second-order AR process, a fact that may be exploited 
later. Various statistical properties of the signal strength process {X[k]} that may be used in 
sequels may also be derived. 
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[0068] A model for hysteresis-based hard handoff mechanisms including a new formulation 
of the handoff problem and the possibility of using a Monte Carlo simulation to obtain 
handoff performance measures will now be discussed. The MS may make a handoff decision 
based on the processed pilot signal strengths from the two candidate base stations. In other 
words, a handoff decision made at time t may be based on the vector X(t) = (Xfl), X/t)) in 
two-dimensional space 1Z 2 . A large class of so-called hysteresis-based handoff algorithms 

maybe characterized in terms of a partition of Tl 1 into three regions (shown in Figures 2, 3, 

and 4): 

I: the assignment region for BS t , 

J: the assignment region for BSj, 

H: the hysteresis region, where the MS may be assigned either to BSi or BSj. 

[0069] A handoff decision at time t > 0 may be made as follows: 

1 . If X(t) e I then a handoff occurs to BS t unless the MS is already assigned to BSi. 

2. If X(t) € J then a handoff occurs to BSj unless the MS is already assigned to BSj. 

3 . If X(t) € H then no handoff occurs. 

At time t = 0 the MS becomes active and is assigned to BSi if X(0) > X/0) and otherwise to 
BSj. In this case, the hysteresis level is effectively zero and there are only two assignment 
regions: 

lo = {(xj; xj : Xi > xj}: the assignment region for BSi at time 0, 
Jo = {(xy, X() : xi < Xj}: the assignment region for BSj at time 0. 
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[0070] With the above formulation of hysteresis-based handoff, some important special cases 
may be considered. A class of handoff algorithms preferably operates on the basis of the 
relative processed signal strength from BSi and BSf. 

X(t)=X i (t)-X j (t) (19) 

[0071] A handoff of the mobile user from base station BSt to base station BSj may occur 
when the relative processed signal strength X(t) falls below the value -hj. Conversely, a 
handoff from base station BSj to base station BSi may occur when X(t) exceeds the value h 
The values h t and hj may be known as the hysteresis levels associated with base stations BSi 
and BSj, respectively. For handoff algorithms based on relative signal strength, the 
assignment regions may be given as follows: 
I = {(*,■; Xi) : xt > xj + k} 

J = {(*,; Xj)_: Xi < xj + hi} 

H = {(x j; x.) : x, - h t < jc, < x, + hi} = K 2 - (J U J) 

These regions are illustrated in Figure 2. 
[0072] Some handoff mechanisms combine relative signal strength with an absolute 

threshold on the required pilot signal strength. In such mechanisms, a handoff from BSi to BSj 
may occur at time t if and only if X(t) < hj and Xft) < Th t , where Th t is the absolute threshold 

on Xi(t). Conversely, a handoff from BSj to BS t may occur at time t if and only if X(t) > hi and 

Xj(t) < Thj. The two-dimensional assignment regions for this handoff algorithm may become: 
I = {(xj; xi) : Xi > xj + hj and Xj<Thj} 
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J = {(*,•; Xii: Xj < xj - hi and x, < Thi) 
H = ft 2 - (I U J) 

[0073] These regions are illustrated in Figure 3. Different specifications of the assignment 
regions lead to different handoff algorithms. A nonlinear assignment region set is illustrated 
in Figure 4. 

[0074] The handoff behavior of an MS traveling along a given trajectory may be viewed as a 
two-dimensional level-crossing problem. Handoff events may occur when the process X(t) 
enters the assignment regions I and J. The level-crossing analogy may be more readily 
visualized in the case of handoff based on relative signal strength. In this case, the handoff 
decision may be based on the one-dimensional relative signal strength X(t). Handoff events 
correspond to level-crossings of the process X(t) with respect to the hysteresis levels h t and - 
hj. Thus, handoffs from BS t to BSj may correspond to downcrossings of -Ay, while handoffs 
from BSj to BSi may correspond to upcrossings of h h as can be seen in Figure 5. 

[0075] The handoff analysis of Vijayan and Holtzman (discussed in the background section) 
is based on the level-crossing analogy. In the asymptotic regime, h t = hj oo, the point 

process formed by the set of upcrossings of h and that formed by the set of downcrossings of 
-hj approach independent Poisson processes. Based on this result, expressions for handoff 
performance metrics such as handoff probability and cell assignment probability may be 
obtained. However, as they point out, the Poisson process assumption is only asymptotically 
true for large h and hj and is not accurate for smaller values of the hysteresis levels. 
[0076] In practice, as the MS moves along a given trajectory, it samples and processes the 
signal strengths from the candidate base stations at discrete time instants. Level crossings 
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between sampling instants do not necessarily lead to actual handoff events. Analysis may rely 
on the assumption that the sampling interval is sufficiently small such that, at most, one level- 
crossing event can occur between sampling instants. Thus, the analysis based on the 
continuous-time model may not capture the effect of varying the sampling interval on handoff 
behavior. 

[0077] In the discrete-time model, the relative processed signal strength between base 
stations BSi and BSj may be given by: 

X[k] = X\k]-Xj[k],k>0. (20) 

Correspondingly, we define 

m[Jfc] = ii % [ik]-« y [*],F[*] = F < [*]-r y [*],*iO. (21) 
Note that process {XfkJ} satisfies the difference equation obtained from (16) if the subscript / 
is removed , i.e., 

X[k] = (a + b)X[k-l]-ahX[k^ (22) 

For handoff algorithms based on relative signal strength, the cell assignment of the mobile 
may be determined by the value of XfkJ, k > 1 . In particular, if the mobile is assigned to BSi 

at time M, then a handoff to BSj occurs at time k if and only if XfkJ < -h h where hj is the 
hysteresis level associated with BSj. Conversely, if the mobile is assigned to BSj at time k - 1, 
then a handoff to BSi occurs at time k if and only if XfkJ > h u where hi is the hysteresis level 
for BS h 

[0078] A Monte Carlo sample path of handoff along a trajectory may be generated 
straightforwardly using the discrete-time model described above. The difference equation 
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characterized by (16) - (18) may determine the processed signal strength {XifkJ} . From (6) 
and (2), the mean pilot signal strength mjfk] may be given as follows: 



The process { Vjfk]} may be simulated simply as a sequence of independent, zero-mean 
Gaussian random variables with standard deviation av,. The processed signal strength, 
{ VifkJ } , may be generated similarly. The relative signal strength process { VfkJ } , given by 
(20), may completely determine the handoff behavior of the mobile along the trajectory. 
[0079] For a given straight-line trajectory, the handoff performance measures of interest 



• Assignment probabilities: Let P t [k] and Pjfk], respectively, denote the probabilities 
that the mobile is assigned to BSi and BSj at time k. 

• Handoff probabilities: Let PyfkJ and Pjifk], respectively, denote the probabilities of a 
handoff from BSi to BSj and a handoff from BSj to BSi at time k. 

• Mean number of handoffs: Let N ho denote the number of handoffs between BSi and 
BSj over the given trajectory. 

• Crossover point: Let C denote the crossover point, which is defined as the point 
along an arbitrary straight-line trajectory at which the probability of the MS being 
assigned to BSi drops below 0:5, given that the probability of assignment to BSi at 
time t = 0 is greater than 0:5. 

[0080] To estimate these performance measures using Monte Carlo simulation, N sample 
paths of the process {XfkJ} may be generated. For the nth sample path, define: 



m\k] 



= v / - j u f log(|a-b,+*v/ I |). 



(23) 



include: 
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[0081] Similarly, Af [k] = 1 if the mobile is assigned to BSj and is equal to zero otherwise. 
Since, for this embodiment, it is assumed that the mobile is assigned to exactly one of the two 
base stations while traversing an arbitrary straight-line trajectory, Af [k] + Af [k] = \. The 
Monte Carlo estimate of the cell assignment probability to BS> at time k may then be given by: 

w „=i 

Under the assumption discussed above, we have P t [k] + Pj [&] = !• 
[0082] Let 4" ) W = 4 W ^- 1 ]4" ) W- Then 4 n) W = 1 ifandonlyifahandofffromitf, 
to BSj occurs at time k in the nth sample path. Similarly, Af [k] = Af [k-\] Af [k] 
indicates a handoff from BSj to BS t at time k in the nth sample path. The Monte Carlo 
estimate of the cell handoff probability from BSt to BSj at time k may be: 

IV n=\ 

[0083] A similar expression may hold for the estimate, P f [k] , of the cell handoff probability 
from BSj to BS t . Let N$ denote the number of cell handoffs that occur over the nth sample 
path of the process {*[*]}. The Monte Carlo estimate of the mean number of handoffs may 
be expressed as: 

*-4t*s?-t(4tWD (27) 

An estimate for the crossover point of the trajectory while the MS travels from cell (/) to cell 
(j) may be given by 

C = argnun{0< it <M :/|[fc]< 0.5} (28) 
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where it is assumed that P. [0] > 0.5 . 

[0084] Analytical expressions for the handoff performance measures in discrete-time and an 
outline for an efficient numerical procedure for computing the handoff performance measures 
will now be discussed. These analytical expressions for the handoff performance measures 
may be derived based on the discrete-time model just discussed. 

[0085] First, assignment probabilities will be considered. The sampled and processed signal 
strength measurements from BSi and BSj may be represented by the vector process X[k] = 
(Xi[k], Xj[k]). Let X[k\r] denote the sequence consisting of the r most recent values of {X[k]} 

up to and including time where 1 < r < k + 1: 

X[k\r] = (X[k-r+l], X[k-r+2], • • • , X[k]) (29) 

[0086] At any time k, the vector X[k] falls in exactly one of the three assignment regions I, J, 
and H. Cell assignment of an MS along a given trajectory may be characterized as follows: 

For k > 1, if the MS is assigned to BSi at time k — 1, a handoff to cell j may occur at time k if 

and only if X[k] G J. Conversely, if the MS is assigned to BSj at time k — 1, a handoff to cell 

i may occur at time k if and only if X[k] G I. Cell assignment may also be characterized as 
follows: 

[0087] Proposition 1: At time k>l 9 the MS is assigned to BSi if and only if X[k\k+1] G I 0 H* 
or X[k\r] G IH r_y for some r; 1 < r < k. Similarly, the MS is assigned to BSj if and only if 
X[k\k+1] G J 0 H* or X[k\r] G JH r " 7 for some r, where 1 < r < k. 
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[0088] Proof: Suppose X[k\r] G IH r / . This condition is equivalent to X[k-r+l] € I and X[/] 
G H, / = k - r + 2, ■ • • , k. Clearly, the MS is assigned to BSi at time k-r + 1. Since no 

handoff occurs at times / = k - r + 2, . . . , k, the MS must also be assigned to BS t at times / = k 
-r + 2,...,k. Similarly, if X[k\ k+1] G I 0 H k , then the MS must be assigned to BS t at time k. 

[0089] Conversely, assume that the MS is assigned to BS t at time k. Then, X[/] G I for some 

1 < / < k or X[/] G lo for / = 0. Let /* be the largest such /. Since the MS is assigned to BS t 

at time k, it must be that X/ G H; m = /*;..., k. Setting r = k-l*+l, we have the desired 

result. The case of assignment of MS to BSj at time k may be proved similarly. 
[0090] Let Ei[k] denote the event that the mobile is assigned to base station BSi at time k. 
Ej[k] may be referred as the cell i assignment event at time k. By virtue of Proposition 1, E[k\ 
may be expressed as follows: 

£ i [it] = {x[A:/A: + l]eI 0 H i }uU{x[A:/r]€lH r - 1 } (30) . 

For convenience, the following notation is introduced: 

p k ( AH'" 1 ) = Pr {X [k I r] e AH'" 1 } , 

where A may be I, l 0 , J, or J 0 . Noting that the events on the right-hand side of (30) may be 
mutually exclusive, we have the following: 
[0091] Proposition 2: The cell assignment probabilities may be expressed as: 

P,[k] = *(*,[*]) = Pk (I 0 H* ) + £ Pk (IH-), (31) 
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Pj [k] = Pr ( Ej [k]) = p k ( J 0 H* ) + ± A K" ) (32) 
10092] Proposition 2 provides a characterization of the cell assignment probabilities under 

the discrete-time handoff model. 
[00931 Handoff probabilities will now be considered. Let Eft] denote the event that at time 
* the cell assignment changes from cell i to cell;, i.e., a handoff may occur from cell i to cell 
j. Similarly, Ej{k] denotes the event that a handoff may occur from cell; to cell i at time k. 
Ite cell handoff events maybe expressed as E M in termS ° f the Ce " 38518111116111 

events defined earlier as follows: 

Ey[k) = E,[k-1] n Ej[k] and E }l {k] = Ej[k-1] n k > 1 (33) 

Analogous to Proposition 2, the conditions for a handoff event to occur may be expressed as 
follows. 

[0094] Proposition 3: For k > 1, a handoff from i to; may occur at time k if and only if X[*] 
6 J and X[k-lV-l] € Ih- for some r; 2 < r < W or X[M\*] € *B» Similarly, a handoff 
from A$ to A* may occur at time k if and only if**] € I and X[fc-lV-l] 6 JH' 2 for some r. 



2<r<fc+lorX[MW|e JoH*" 1 . 



[00951 Hence, the handoff events and £tf *] may be expressed as follows: 

^W = {x[H/c + l]eI 0 H-j}uy{x[^r]eIH-j} (33) 

£ J^ = {x[a^l]€J„H-l}uU{x[Hr] e JH-l},^l (34) 
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Since the events on the right-hand sides of (34) and (35) are mutually exclusive, the following 
characterization of the cell handoff probabilities, analogous to Proposition 2, may be made: 
10096] Proposition 4: The cell handoff probabilities maybe expressed as follows: 

^W = Pr(^W) = p t (l 0 H*--j) + tp t (lH-j)k 06) 



P,W = Pr(£,l^]) = ^(j„H t - 1 l) + t/'*K 2 l) s **1 (37) 

[0097] An important performance measure to estimate the context of a cell assignment model 
is the mean number of handoffs that occur as the mobile moves along a given trajectory. Let 
N h0 denote the number of handoffs that occur for a mobile moving along a given trajectory. 
Let 1, denote the indicator function for the event E. The number of handoffs may then be 
expressed as follows: 

(38) 

[0098] Hence, the mean number of handoffs may be given by 

£["J=£H E <M) +Pr ( E *M)) <39) 

k=\ 

[0099] Recursive procedure for handoff analysis will now be considered. For handoff 
algorithms based on relative signal strength, the handoff behavior of the MS may map to the 
one-dimensional process X[*]. In this case, the assignment regions may be intervals of the 
real-line: / = [h& oo],. J= (-oo; h,l and H- (-hf, V). The assignment regions at system 
initialization time may be I 0 = (0,oo) and J 0 = (-oo; 0). The expressions derived earlier for 
assignment and handoff probability map straightforwardly to the one-dimensional case. In the 
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one-dimensional case, equation (31) may express the cell assignment probability to BSi at 
time k as a summation of probabilities of the form pk{ttf A ). 
[00100] A recursive procedure for evaluating these probabilities will now be developed. For 

convenience, X[k] may be denoted by X*. A sequence of bivariate functions {gkixk-itXk), k > 

1 } may be defined as follows: 

gi(*o>X\) =A>Xi(*o>*i) (40) 
Sii^i) = [ gii^^f^ix^x^x^dx, (41) 

*1q 

Si { X k-\> X k) = [ fak-2 X k-\ { X k-2> X k-\) Akl X *-2 X *-l { X k I X Jt-2' X ifc-l)^-2 + 

f / , x < 42 > 

= J„S*-1 (**-2>**-l)Akl X *-2 X *-1 (** I **-2.**-l)<&*-2. * * 3 

[00101] Here, ./x^x* is the joint density of {X [k- l],X[ A:]} and 

Akl x *-2 x *-i (■** I x k-2> x k-\ ) i s me conditional density of X[k], given X[k-l] and X[&-2]. 
Expressions for the conditional joint densities may be derived by those skilled in the art. The 
notation £ denotes the standard Riemann integral over the set A. The bivariate functions g* 
have the following useful property: 
[00102] Lemma 1: Fork > 2, 

G k = I \ H g k (x k _ x ,x k )dx k dx k _ x =P k {l 0 H k ) + tp k (lH r - 1 ) (43) 

[00103] With these definitions, the probability of cell assignment to BSi may be expressed as 

follows (the corresponding result for the probability of assignment to BSj being analogous): 
[00104] Proposition 5: 

p i M = Pi i 1 ) + M IH )+H S* (**-. . x * >M**-. . * * 2 > (44) 
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/J[i] = aM+a(VO, < 45 > 

fl[0] = /%(/.). (46) 
[00105] The probabilities />*(I) and /?*(IH) may be given as follows: 

a(/)={AM* 4 ,**0 (47) 

^* ( 7i/ ) = I JL ■ /k *- iX * (**-■ 5 x " K^*-' ' * - 1 • (48) 



[00106] Proof: Eqs. (46) and (45) follow directly from Proposition 2 when k = 0 and k = 1, 
respectively. For k > 2, Eq. (44) follows from Proposition 2 and Lemma 1. 

[00107] The handoff probability from BSi and BSj may be expressed similarly in terms of the 
functions gk (the corresponding result for the probability of handoff from BSj to BSi being 
analogous): 

[00108] Proposition 6: 

Pij M = A W £ [gk (**-,.** >M**-„ * * 2 ' ( 49 > 

A-[1] = a(V)- (50) 
[00109J The probability piiJJ) may be given as follows: 

Ptty]= JjAAfw)*^.^ 1 ' ( 51 ) 
[001 10] A methodology for characterizing handoff performance over arbitrary mobile 
trajectories in a cellular network having a regular hexagonal cell pattern will now be 
disclosed. It is envisioned that this approach may be extended to networks with other types of 
regular cell structure. Figure 6 illustrates an example of a mobile trajectory inside a wireless 
network covered by hexagonal cells. In the hard handoff scheme, at a given location, only the 
two base stations closest to the mobile unit are involved in the cell assignment process. (Of 
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course, in a soft handoff scheme, more that two base stations would be involved). The two 
base stations preferably define a rhombus in which the major diagonal is a line segment 
joining the two base stations. This rhombus may be referred to as a minimum geometric 
structure (MGS). 

[00111] For the purposes of modeling and analyzing hard handoff performance, an arbitrary 
trajectory may be mapped into a single instance of an MGS. In order to accomplish this, the 
concept of a reduced geometric structure (RGS) is introduced first. As shown in Figure 7, an 
RGS may be a star-shaped object consisting of a hexagonal cell at the center and six 
triangular regions subtending the edges of the hexagon. The mapping of the cellular network 
into a single MGS may be performed in two steps: 

1 . translation of all path segments into a single reduced geometric structure (RGS); 

and 

2. rotation of all path segments within the RGS into a single MGS. 

[001 12[ In the hexagonal cellular network, the RGS preferably has the shape of a star, with a 
base station at the center and six neighboring base stations situated at the vertices of the star. 
The RGS may be viewed as consisting of six MGSs, one corresponding to each of the 
vertices. 

[001 13] The mapping of a mobile trajectory in a cellular network into an RGS is illustrated in 
Figures 6 and 7. hi Figure 6, base station BS, lies at the center of an RGS that is chosen as the 
base RGS. The mobile trajectory begins at location r 0 in the RGS centered at BS\ and enters 
the RGS centered at BSi, at location r 2 . Since the path segment r 0 r 2 lies entirely within the 
base RGS, no translation may be required. The path segment r 2 r 5 lies within the RGS 
centered at BS 2 and should be translated to the base RGS. This may be done by translating the 
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RGS centered at BS3 to the RGS centered at BSj. Figure 7 shows the segments of the original 
trajectory that have been translated to the base RGS. For example, the path segment r 2 r 5 , 
originally lying within the RGS centered at BSi in Figure 6 has been translated to its 
corresponding position within the base RGS in Figure 7. Similarly, the path segment r 5 r 6y 
originally lying within the RGS centered at BS4 has been translated to its corresponding 
position in the base RGS in Figure 7. Thus, Figure 7 shows a representation of the original 
mobile trajectory in the cellular network that has been mapped to a single base RGS via 
translation. 

[00114] The trajectory segments shown in the base RGS of Figure 7 may be further mapped 
into an MGS via rotations. The RGS in Figure 7 consists of six rhombi corresponding to the 
vertices of the RGS. The rhombus corresponding to the leftmost vertex of the RGS to serve 
as the base MGS may be chosen. The remaining five MGSs may each be rotated in turn to 
coincide with the base MGS. Within an MGS, a trajectory that passes from one boundary to 
another may be approximated by a line segment. As shown in Figure 9, the segment r\ri may 
be characterized uniquely by the crossing point located at the distance d from Oq and the 
angle a, which it forms with (O0O2). A set of line segments within the MGS may be 

generated by the parameter pair {d\ a) over the range (0; </ max ) x (0; n) using a particular step 

pair (8 d ; 5 a ). The overall handoff performance within a cellular network may be 
characterized in terms of the mean number of handoffs along each segment in the set of 
generated line segments. 
[00115] Numerical evaluations of handoff performance obtained using embodiments of the 
present invention will now be discussed. Comparisons were made with results obtained from 
Monte Carlo simulation as well as the alternative analytical approaches mentioned earlier. 
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First we evaluate the cell assignment and handoff probabilities along two representative 
trajectories Tr\ and Tr 2 shown in Figure 9. Within the minimum geometric structure of 
Figure 9, the trajectories Tr\ and Tr 2 are characterized by the parameter pairs (di; en) = (0; 0) 
and (d 2 ; a 2 ) = (1000; n/4), respectively. Also evaluated is the impact of the hysteresis 
parameter h on the mean number of handoffs and the crossover point for the two trajectories. 
Finally, handoff performance over the set of straight-line trajectories generated within a 
minimum geometric structure is evaluated. 

[001 16] Numerical results may be obtained using an embodiment implemented in MATLAB 
of the recursive procedure described earlier. The implementation represents the bivariate 
function gk(*k-ir*k) as a matrix over a finite grid of points in two-dimensional space. The 
integrals that appear in the recursive procedure are approximated by trapezoidal integration. 
The numerical procedure may be made arbitrarily more accurate by using a larger matrix 
corresponding to a grid of higher granularity, but at the expense of greater storage 
requirement and higher computational cost. The procedure computes assignment probabilities 
to an error of less than 10" 6 . The overall computational complexity of the procedure at each 
iteration is 0(h 2 ), where h is the hysteresis parameter. It is believed that this implementation 
is more efficient than Monte Carlo simulation for an equivalent level of accuracy. 

[00117] For the illustrative numerical computations presented, the main system parameters 
were set as follows: D = 2000 m, vj = 0 dB, u { = 30 dB, aw { = 6 dB, d 0 = 20 m, d s = 1 m. The 
decay distance of the smoothing window d zv was set to 10 m. For comparison purposes, 
numerical results obtained using the following alternative approaches to the handoff 
performance evaluation are also shown: 
• Monte Carlo simulation; 
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• Poisson approximation model of Vijayan and Holtzman; and 

• Approximate discrete-time model of Zhang and Holtzman. 

[00118] Figures 10 shows the probability that the mobile station is assigned to base station BSi 
as a function of the sampled distance along the trajectory Tr, when the hysteresis parameter is 
set to h = 1 . The curve obtained from the present model matches closely that obtained from 
the Monte Carlo simulation. Since the mobile begins at BS t , the assignment probability is 
initially 1 and then decreases monotonically to 0 as it approaches BSj along the straight-line 
trajectory. The cell assignment probability along Tr, obtained using the Vijayan-Holtzman 
model is also plotted in Figure 10. The assignment probabilities obtained from the Vijayan- 
Holtzman model show inaccuracies, and a bias may be seen in the crossover point. For 
relatively small hysteresis levels, the Vijayan-Holtzman model does not appear to provide a 
good estimate for the crossover point and the mean number of handoffs. Similar observations 
may be made with respect to Figure 11, which shows the assignment probability to base 
station BS< as the mobile station moves along the trajectory lr 2 for a hysteresis value of h = 1 . 
[00119] Observe in Figure 10 that the crossover point is approximately equal to the halfway 
point along the trajectory, i.e., 1000 m. The ideal assignment probability curve along the 
trajectory Tr, is a step function with a step discontinuity from the level 1 to the level 0 at the 
halfway point. As the hysteresis level h is increased, the assignment probability curve more 
closely approximates a step function as may be seen in Figure 12. However, in this case, as 
the value of the crossover point increases, a larger handoff delay is implied. Therein lies a 
critical tradeoff in dimensioning the parameter h for hysteresis-based handoff algorithms. 
[00120] Comparing the Vijayan-Holtzman results in Figures 1 0 and 1 2, the results may 
become more accurate as h is increased. As noted by Vijayan and Holtzman, the Poisson 
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model is asymptotically accurate as the value of the parameter h is increased. However, as 
shown in Figure 13, even when h = 10, a discrepancy may be observed between the 
assignment probability curve from the Vijayan-Holtzman model and that from simulation and 
this analysis. Also note that there is significant error in the first couple of hundred meters, 
which is approximately the radius of a microcell. This is most likely because the model does 
not take into account the transient behavior due to system initialization when the MS becomes 
active. This effect is more pronounced for larger h. 

[00121] Figure 14 shows the probability of the mobile station being handed off from BSi to 
BSj as it moves along trajectory Tr\ when the hysteresis parameter is set to h = 3. The 
simulation results are shown along with 95% confidence intervals. The curve obtained from 
our analysis agrees closely with the simulation curve, whereas the curve obtained from the 
Vijayan-Holtzman model underestimates the handoff probabilities. 

[00122] Similar observations may be made with respect to Figure 1 5, which shows the 

handoff probabilities for the mobile station moving along trajectory Tr 2 . For both trajectories, 
the handoff probability increases to a maximum value near the crossover point. It should 
further be noted that for small values of the handoff probabilities, a large number of 
simulations may be required, making the present analytical approach an attractive alternative 
to simulation. 

[00123] A main performance measures used in the design and optimization of the 2G wireless 
networks are the mean number of handoffs and the crossover point. In Figure 17, the mean 
number of handoffs for a mobile station moving along trajectory lr\ is plotted as a function of 
the hysteresis parameter h using the four methods discussed above. The present analytical 
results agree closely with those obtained from Monte Carlo simulation. Note that the Zhang- 
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Holtzman model overestimates the mean number of handoffs over almost the entire range of 
hysteresis values. The Vijayan-Holtzman model underestimates the mean number of handoffs 
when h is less than 5, while for h>5, this model shows close agreement with the present 
model and Monte Carlo simulation. However, as shown in Figure 16, the Vijayan-Holtzman 
model does not capture the transient handoff behavior at the beginning of the MS trajectory 
for the first couple of hundred meters. 
[001 24] Figure 1 9 plots the crossover point as a function of the hysteresis parameter for the 
mobile station moving along trajectory Tn. Figure 20 shows the corresponding results for 
Tr 2 . In both figures, results from simulation and the analytical model are shown. The general 
trend is that the crossover point increases as the hysteresis level increases. Thus, larger 
hysteresis levels result in larger delays in handoff. On the other hand, as observed earlier, the 
mean number of handoffs decreases with increasing hysteresis level. Depending on the 
application, different tradeoffs between the mean number of handoffs and the crossover point 
may be desirable. The hysteresis parameter should be chosen accordingly. 
[00125] In general, it is desirable to minimize the number of handoffs while having as small a 
handoff delay as possible. A small handoff delay may be particularly important for the design 
and dimensioning of 3G systems, since the allocated bandwidth is significantly larger than for 
2G systems. As seen earlier, small handoff delay may be achieved by using small values for 
the hysteresis parameter. On the other hand, small hysteresis levels may imply a larger 
number of handoffs. The tradeoff between crossover point and mean number of handoffs for 
trajectories Tn and Tr 2 is illustrated explicitly in Figures 21 and 22, respectively. In both 
figures, the optimum operating point will typically lie in the neighborhood of the knee of the 
curve. 
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[00126] Using the minimum geometric structure introduced earlier (see Figure 9), the mean 
number of handoffs for a cellular network with hexagonal structure may be characterized as a 
surface. In particular, we characterize the mean number of handoffs over a set of line 
segments within a minimum geometric structure {QoOxOiOi). 

[00127] The set of line segments for the set of line segments generated by the following set of 
relations may be shown as: 

{(d,«)6(0,2000)x(0,*);(rf e ,a) = *(^,5 a ),* = ai,-} (53) 

[001 28] Figures 23 and 24 show the mean number of handoffs obtained from this handoff 
model over the set of straight-line trajectories within a minimum geometric structure when the 
hysteresis parameter h is set to 1 and 3, respectively. For these figures, the set of segments 

was generated using (53) with ,S a ) = (lOO,^) • In Figures 23 and 24, the mean number 

of handoffs may achieve a maximum when d c = 1000 and a = n2, which corresponds to the 
trajectory 7V 3 along the cell boundary between 55, and BSj. For this particular trajectory, the 
mean number of handoffs is about 12 when h = 1 and about 6 when h = 3. As h is increased, 
the entire surface becomes lower, though not in a uniform manner. 
[00129] Each surface summarizes the handoff performance over the cellular network for a 
particular hysteresis value. As discussed earlier, an arbitrary trajectory in the network may be 
obtained by approximating the trajectory with a set of straight-line trajectories within a 
minimum geometric structure. The mean number of handoffs along the trajectory may then 
be obtained using surfaces such as those shown in Figures 23 and 24. These surfaces provide 
a convenient means for visualizing the overall signaling load due to handoffs in the network. 
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[00130] A main contribution of the present invention is a rigorous discrete-time framework 
for modeling and analyzing handoff performance in cellular networks. Based on the discrete- 
time framework, expressions may be derived for the probabilities of cell assignment and 
handoff along arbitrary straight-line trajectories. An efficient numerical procedure may then 
be developed for accurately evaluating the performance measures for handoff algorithms 
based on relative signal strength. By introducing the concept of a minimum geometric 
structure, a characterization of handoff performance over a wide range of possible mobile 
trajectories is developed. This invention finds direct application to the design and 
dimensioning of handoff algorithms in general and more specifically in hard handoff 
algorithms in cellular networks based on the GSM and GPRS standards. 
[00131] In recent years, some novel handoff algorithms have been proposed in the research 
literature. For example, some algorithms formulate a locally optimal handoff algorithm to 
obtain the best tradeoff between expected number of service failures and expected number of 
handoffs. Other algorithms suggest the use of neural networks for handoff algorithms. Such 
algorithms may result in nonlinear assignment regions (cf. Figure 4). In principle, nonlinear 
assignment regions are accommodated in the presently disclosed discrete-time framework, but 
the numerical evaluation of the assignment and handoff probabilities in this case may be 
computationally challenging. The same remark applies to handoff algorithms that use both 
relative and absolute signal strengths (cf. Figure 3). 
[00132] The numerical procedure developed for the present handoff mechanisms may be 
efficient for relatively small values of the hysteresis parameter h, but its complexity grows as 
0(h 2 ). This complexity is due to the second-order Markov property of the processed signal 
strength {A[Jfc]}. By contrast, the raw signal strength {Y[k]} maybe characterized as a first- 
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order Markov chain. The higher order of the Markov chain {X[k]} may be due to the 
averaging filter/ av [A;], which is used to eliminate fast fading. An alternative technique may be 
to eliminate fast fading while preserving the first-order Markov property of the raw signal 
strength. In this case, the handoff analysis may be performed more efficiently than with 
conventional averaging. 

[00133] Although the present invention has been described illustratively in the case of hard 
handoff, the discrete-time handoff model may be extended to the case of soft handoff. For 
example, it is envisioned that the present invention may be applied to soft handoffs and the 
design of packet-switched CDMA networks that exploit macro-diversity. 

[00134] Another illustrative embodiment of the invention which focuses specifically on soft 
handoff in IS-95 will be presented. This software handoff mechanism may be characterized 
by three basic parameters: an add threshold T a , a drop threshold T d , and a drop timer r d . As 
the mobile unit moves in the network, it may sample received pilot signal strengths from base 
stations in its close vicinity, which may form a neighbor set, N k , where k is the discrete-time 
sampling index. Among the base stations in N k , some belong to the active set, A k , and 
others to the candidate set, C k . The candidate set is the set of neighbor base stations not in the 
active set with pilot signals of sufficient strength for the mobile unit to demodulate the 
associated forward traffic channels. A base station in the candidate set enters the active set if 
its pilot signal strength exceeds T a . A base station currently in the active set leaves the active 
set if its pilot signal strength drops below T d for a period of time exceeding x d . 

[00135] A cell assignment model will now be discussed. The assignment of cell i to the 
mobile unit f s active set may be determined by the time-evolution of the received signal 
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strength process Y.(t). In IS-95 CDMA networks, cell i is preferably added to the mobile's 
active set when the value of Y.(t) exceeds the add threshold T a . Cell i may then be removed 
from the active set when Y.(t) falls below the drop threshold T d for an interval of at least r d 
seconds. In the discrete-time framework, the cell assignment decision may be based on the 
process Y^k]. Cell i may be dropped from the active set when the value of Y^k] is less than 
T d for M consecutive discrete-time instants, where M = \r d lr s \ The parameters T a , T d , 
and M may determine the soft handoff algorithm. In general, these parameters may also 
depend on the cell i . 

[00136] An illustrative soft handoff embodiment of the present invention will now be 
discussed. A real line may be partitioned into three disjoint intervals defined as follows: 

A = [T a ^\ H = (T d ,T a ), D = {-n,T d \ (54) 
where A is the assignment region, H is the hysteresis region, and D is the drop region. If 
Y^e A , then cell i should be assigned to the active set at time k , i.e., i e A k . Let Y.[k \ r] 
denote the vector consisting of the r most recent values of Y t \k\ up to and including time k : 
I5[*\r]=fc[*-r + ll... i i;M (55) 

[00137] A soft handoff rule at time k > 1 may be expressed as follows: If cell i is not in the 
active set at time k - 1 , it may enter the active set at time k if and only if Y t [k] e A . 
Conversely, if cell i is not in the active set at time k - 1 , it may leave the active set at time k 
if and only if Y t [k \ M] g D m . Cell / may be assumed not to be in the active set at time k = 0. 

[00138] It is convenient to characterize the received signal strength in terms of the regions A , 
H , and D . Specifically, the received signal Y^k] may be mapped to a letter S^k] defined on 
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the alphabet S = {A,H,D}. Let S" denote the set of strings of length n over the alphabet S . 
The cardinality of S" is preferably 3" . Let S + denote the set of all finite length strings over 
the alphabet S , i.e., S + = <J* =l S" . 
[00139] Let Ei[k] denote the event that cell i is assigned to the active set at time k . The 
complementary event, Ei[k] should be the event that at time k cell i is not assigned to the 
active set. Then the cell non-assignment event, Ei[k], preferably occurs if and only if one of 
the following conditions holds: 

Y[k]iA for n = l k, or (56) 

Y.[n \M] e D M for some n, M < n < k and Yfr] £ A for n < 1 < k. (57) 
[00140] Let B k denote the set of strings of length k over the alphabet S that do not contain 
the character A . B k may be expressed as B k = {H,Df . The first condition for cell non- 
assignment may be formulated in terms of the string S t [k \ k] as follows: 

S,[k\k]eB k (58) 
[00141] Let H k denote the set of strings of length k that do not contain the sub-string D M or 
the character A . The set H k may be expressed as follows: 

H k =B k Kj {h,HD,HD\...,HD m - x \ (59) 
The second condition for cell non-assignment may then be formulated as follows: 
S i [k\r]eD M H r - M for some r, M < r < k. (60) 

The conditions (58) and (59) may be combined to express E> [k] as a union of disjoint events: 
£,[*] = \ *] € B k } u (J \s,[k \ r] € D M H r _ u } (61) 
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For a string s e S + of length \s\<k, the following notation may be introduced: 

[00142] Thus, p kJ (s) may represent the probability that the last|s| values of the sequence 
Y\k \ k] fall in the regions defined by the strings . The probability that cell i is assigned to 
the active set may be denoted by P\k] = P{E i [k]} . Similarly, the probability that cell i is not 
assigned to the active set may be denoted by Pi[k] = Now, 
/>[*]+ P, [fc] = l for all* >0. 

[00143] The characterization of cell non-assignment in (61) may imply the following result. 

[00144] Proposition 7: 

aw-i^w+i Za» (63) 

[00145] Proposition 7 provides the basis for the numerical evaluation of the cell assignment 
probabilities. From the cell assignment probabilities, a number of handoff performance 
measures may be obtained. The outage probability is the probability that the active set is 
empty and may be given by: 

p.M=rF*i (64) 

ieN v 

where N " is the set of neighbor cells. The mean size of the active set or macrodiversity 
degree may be given by: 

N*=4iE*i (65) 

A k=li€N k 
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[00146] An efficient recursive procedure for evaluating the non-assignment probabilities 
given in Proposition 7 will now be presented. To that end, a family of M - 1 sequences of 
functions may be defined as: 

h\ m] :D^9{\ k>\,\<m<M-l, (66) 

where ^ may denote the real line, as follows: 

hiky)=lf k _ x (z)f k {y\z)dz, 

m = 2,...,M-1, 

with the initialization condition: 

hl m] (y) =0 for all y e D if m < k. (68) 
[00147] Here, f k (x) may denote the marginal density of Y[k] and f k {x\z) may denote the 
conditional density of Y[k] given that Y[k - 1] = z . These densities are preferably Gaussian. 
By further observing that Y[k] is most likely a first-order Markov chain, one may show that: 

pXd m )= ih [ r ] iy)dy- (69) 

[00148] Next, we define a sequence of functions: 

g k : H -> 9T and g l ™ ] : D -> 9T , (70) 

for m = 1,. . .,M - 1 . For k > 2 , qr t : Z) -» 9T may be defined as follows: 
Next, for > 1 define: 

W= i^., W/*(y I *)* + sl 1] M (72) 

where 
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(73) 



and 



(74) 



for r = 2, . . .,M - 1 . The initialization conditions may be: 



*.M = /• W«»d st" 1 M = 0 for k < m. 



(75) 



With the above definitions, the cell non-assignment probabilities may be expressed as 
follows. 

[00149] Proposition 8: 



right-hand-side of (76) is equivalent to the expression given in Proposition 7. In doing so, one 
uses the fact, noted earlier, that Y[k] may be a first-order Markov chain. 

[00151] Figure 25 illustrates a possible cellular network and a mobile trajectory along the 
straight line directed from BS l to BS % . In this example, the distance between two 
neighboring base stations is D = 2000 m . The candidate set is C = {1,2,3,4} . Cell assignment 
probabilities and the outage probabilities as functions of the discrete-time parameter k for a 
mobile unit traveling along trajectory in Figure 1 may be calculated using the present 
invention. These analysis may be performed using an implementation of the recursive 
procedure of Proposition 8 in MATLAB. 

[00152] In order to study the mean macrodiversity degree for arbitrary trajectories in the 
cellular network, a family of straight-line segments directed from a fixed origin O may be 
considered as shown in illustrative Figure 26. In this illustrative example, the segments all 



[00150] Proof: (Sketch) The proof is by induction on k . One shows that for each k , the 
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have length d r . The starting point and orientation of each segment is specified by the pair 
(d,a), where d is the radial distance from the origin 0 and a is the angle of the segment 
with respect to the reference axis OX . Equivalently, the pair (d,a) may be expressed in 
Cartesian coordinates as (d x ,d Y )> wnere d x =d cosa and d Y = d sma . 

[001 53] This illustration of a software handoff embodiment of the present invention shows a 
new discrete-time modeling and analysis technique to evaluate the performance of soft 
handoff analysis for mobile units traveling along an arbitrary trajectory in a cellular CDMA. 
An advantage of this embodiment is that it provides for a recursive computational procedure 
to evaluate the cell assignment probabilities efficiently and accurately. This analytical 
approach is advantageous to computer simulation both in accuracy and in efficiency, 
particularly in calculating outage probabilities. This invention may also be used to dimension 
and optimize soft handoff algorithms for CDMA networks. Numerical results for a cellular 
network with hexagonal cells show that the value of the drop timer parameter, r d , may have a 
significant effect on the mean macrodiversity degree. It is envisioned that this soft handoff 
analysis system could be used to evaluate the achievable macrodiversity gain of cellular 
network architectures for a variety of indoor and outdoor propagation models. 

[00154] Figure 27 shows a flow diagram of an embodiment of an aspect of the present 
invention. This flow diagram represents basic steps that may be taken to practice an 
embodiment of the present invention. One skilled in the art will recognize that other 
sequences may be performed in accordance with the principals taught in this disclosure. At 
step S2710, a description on a cellular network may be created. This description should 
represent the network where handoff decision mechanisms are to be studied. A description as 
per this invention may include elements such as a cellular network geometry; a signal 
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propagation description; a description of network base stations, and a description of network 
mobile units. The network base stations may be further described using elements such as the 
base stations location within said cellular network geometry and the base stations pilot signal 
strength. The mobile unit(s) preferably are capable of moving along a trajectory and sampling 
pilot signal strengths. 

[00155] At step S2720, a reduced geometric structures may be created from the reference 
cellular network description. A representation of a general trajectory for a mobile unit 
traveling in the cellular network may be generated at step S2730. Next, at step S2740, a 
piecwise path may be approximated from the general trajectory representation. Next, a 
handoff mechanism maybe selected at step S2750 and handoff parameters selected at step 
S2760. Examples of handoff mechanisms include hard handoff mechanisms and soft handoff 
mechanisms. Examples of handoff parameters include: hysteresis values; drop timer values; 
averaging parameters; and handoff decision intervals. 

[00156] A discrete-time formulation may then be created at step 2770. This discrete-time 
formulation preferably characterizes handoff behaviors. Finally, handoff performance 
metric(s) may be calculated at step S2780. These metrics may include: an assignment 
probabilities; a handoff probabilities; mean number of handoff value(s); crossover point(s); 
outage probabilities; macrodiversity degrees; and overall signaling loads incurred by the 
handoff mechanism. 

[00157] Figure 28 shows a block diagram of an embodiment of an aspect of the present 
invention. This block diagram represents basic elements that may be used to practice an 
embodiment of the present invention. One skilled in the art will recognize that other elements 
may be used in accordance with the principals taught in this disclosure. The cellular network 
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modeler 2800 preferably includes: a cellular network generator 2810; a reduced geometric 
structure creator 2820; a general trajectory generator 2830; a piecewise path creator 2840; a 
handoff mechanism selector 2850; a handoff parameter selector 2860; a discrete-time 
formulation creator 2870; and a handoff metric calculator 2880. The cellular network 
generator 2810 should be capable of generating a description of a reference cellular network 
to be studied as previously discussed. The reduced geometric structure creator 2820 is 
preferably capable of creating reduced geometric structures from a reference cellular network 
description. The general trajectory generator 2830 is preferably capable of representing a 
general trajectory for a mobile unit in a reference cellular network description. The piecewise 
path creator 2840 is preferably capable of approximating a general trajectory with a piecewise 
path having a multitude of handoff decision points. The handoff mechanism selector 2850 is 
preferably capable of selecting for analysis a handoff mechanism such as a hard handoff 
mechanism or a soft handoff mechanism. The handoff parameter selector 2860 is preferably 
capable of selecting for use in handoff analysis parameters such as hysteresis values; drop 
timer values; averaging parameters; and handoff decision intervals. The discrete-time 
formulation creator 2870 is preferably capable of creating a discrete-time formulation such as 
previously described. And finally, the handoff metric calculator 2880 is preferably capable of 
calculating performance handoff metrics such as: assignment probabilities; handoff 
probabilities; mean number of handoff value(s); crossover point(s); outage probabilities; 
macrodiversity degrees; and overall signaling loads incurred by the handoff mechanism. 
[00158] It is envisioned that each of these elements may be implemented using digital 
hardware with or without microprocessor assistance. In one embodiment, the elements may 
all be modules operating within a single computer. However, to speed up results, the modules 
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may be split into separate computers that are networked. For certain modules, it will likely be 
advantageous to implement their functions using signal processing. 
[00159] The foregoing descriptions of the preferred embodiments of the present invention 
have been presented for purposes of illustration and description. They are not intended to be 
exhaustive or to limit the invention to the precise forms disclosed, and obviously many 
modifications and variations are possible in light of the above teaching. The illustrated 
embodiments were chosen and described in order to best explain the principles of the 
invention and its practical application to thereby enable others skilled in the art to best utilize 
the invention in various embodiments and with various modifications as are suited to the 
particular use contemplated. 
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